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Magnetization of the lunar crust 
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[i] Magnetic fields measured by the satellite Lunar Prospector show large scale features 
resulting from remanently magnetized crust. Vector data synthesized at satellite altitude 
from a spherical harmonic model of the lunar crustal field, and the radial component of the 
magnetometer data, have been used to produce spatially continuous global magnetization 
models for the lunar crust. The magnetization is expressed in terms of localized basis 
functions, with a magnetization solution selected having the smallest root-mean square 
magnetization for a given fit to the data, controlled by a damping parameter. Suites of 
magnetization models for layers with thicknesses between 10 and 50 km are able to 
reproduce much of the input data, with global misfits of less than 0.5 nT (within the 
uncertainties of the data), and some surface field estimates. The magnetization distributions 
show robust magnitudes for a range of model thicknesses and damping parameters, 
however the magnetization direction is unconstrained. These global models suggest that 
magnetized sources of the lunar crust can be represented by a 30 km thick magnetized 
layer. Average magnetization values in magnetized regions are 30-40 mA/m, similar to the 
measured magnetizations of the Apollo samples and significantly weaker than crustal 
magnetizations for Mars and the Earth. These are the first global magnetization models 
for the Moon, providing lower bounds on the magnitude of lunar crustal magnetization in 
the absence of multiple sample returns, and can be used to predict the crustal contribution 
to the lunar magnetic field at a particular location. 

Citation: Carley, R. A., K. A. Whaler, M. E. Purucker, and J. S. Halekas (2012), Magnetization of the lunar crust, J. Geophys. 
Res., 117, E08001, doi: 10.1 029/20 1 1JE003944. 


1. Introduction 

[ 2 ] Since the Apollo missions it has been known that 
while the Moon possesses no global magnetic field, it has 
regions of weak magnetic field resulting from areas of the 
crust that are remanently magnetized. The NASA orbiting 
satellite Lunar Prospector [Binder, 1998] made measure- 
ments of the lunar magnetic field with a magnetometer 
(MAG), and made estimates of the surface field magnitude 
with an Electron Reflectometer (ER), providing the first 
global maps of the lunar magnetic field [Richmond and 
Hood, 2008; Purucker, 2008; Mitchell et al, 2008], 

[ 3 ] The main features of the lunar crustal field are their 
association with the antipodes of the youngest and similarly 
aged large impact basins [Lin etal., 1988; Mitchell et al., 2008; 
Halekas et al., 2001; Hood et al., 2001], and the presence of 
strong isolated anomalies. While the antipodes of the young 
impact basins show very strong magnetic fields, the large 
mare-filled impact basins themselves, and other lunar 
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craters, show magnetic lows compared with their sur- 
roundings [Halekas et al., 2002], with the exception of 
some older (pre-Nectarian/Nectarian), larger (>50 km 
diameter) impact basins which show weak central magnetic 
fields [Halekas et al., 2003; Purucker et al., 2008; Hood, 
2011]. Many of the strong anomalies, such as the Reiner 
Gamma anomaly, are associated with regions of high 
albedo [Richmond et al., 2003], often swirl-like in nature 
[e.g., Hood and Schubert, 1980], suggested to result from 
differential weathering (optical maturation) of the lunar reg- 
olith due to deflection of the solar wind particles by the 
magnetic anomaly. 

[ 4 ] The large scale (>500 km) magnetic anomalies mea- 
sured by satellites imply the presence of relatively strong 
and coherent sources of remanent magnetization. However, 
surface magnetometers showed crustal fields with great 
variety in strength and direction over distances of only a few 
km, suggesting they are associated with localized sources of 
magnetization close to the surface [Dyal et al, 1970]. Large 
natural remanent magnetizations of the returned Apollo sam- 
ples [Runcorn et al, 1970] showed that the main carriers for 
this remanent magnetization were metallic iron particles 
[Fuller, 1974]. Iron can be derived from the bedrock or can be 
added to lunar soils and breccias during impacts [e.g., Housley 
et al, 1973], and the strongest fields are in regions linked 
with impact ejecta [Strangway et al., 1973; Dyal et al, 1974], 

[ 5 ] Interpretation of lunar magnetism relies on the deter- 
mination of the distribution, magnitude and direction of the 
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sources of magnetization. In the absence of extensive 
sampling, or information about possible source locations 
and dimensions, this is made possible with magnetization 
models produced from the observed magnetic field data. 
Inverse models that do not require assumptions about the 
magnetization direction are particularly important for the 
lunar case where directions cannot be assumed. 

[6] The origins of the magnetic fields that resulted in this 
remanent magnetization, and the mechanisms by which the 
remanent magnetization was acquired, are still under dis- 
cussion. The main contenders for the lunar paleofields are an 
ancient lunar dynamo, and processes occurring on the Moon 
as a result of impacts [e.g., Collinson, 1993], with magne- 
tization features resulting from excavation, heating and 
cooling in the presence of a field, or shock magnetization 
(or demagnetization). Given the highly cratered surface of 
the Moon, and the association of many of the lunar mag- 
netic anomalies with impacts or impact ejecta, it is likely 
that impact processes are an important feature of lunar 
magnetism. 

[ 7 ] The strong remanent magnetization measured in the 
Apollo samples [e.g., Fuller and Cisowski, 1987; Garrick- 
Bethell et al., 2009] and the central basin anomalies might 
suggest the existence of a past lunar dynamo [e.g., Hood, 
2011]. While several mechanisms have been proposed for 
dynamo generation [e.g., Stegman etal., 2003; Le Bars et al., 
2011; Dwyer et al, 2011], magnetization in a dynamo 
field may not be able to explain all of the features of the 
lunar crustal fields. The antipodal relationship between 
impact basins and magnetic anomalies has been explained 
by the focusing of impact generated plasmas and amplifi- 
cation of an ambient magnetic field at the antipode [Hood 
and Vickery, 1984; Hood, 1987; Hood and Huang, 1991; 
Hood and Artemieva, 2008], with magnetization occurring 
through low pressure shocks [Gattacceca et al., 2010]. 
Hypervelocity impacts are also capable of producing large 
transient magnetic fields [Smka, 1977; Crawford and Schultz, 
1988], 

[s] Whatever the magnetic fields responsible for produc- 
ing the remanent crustal fields, they must be able to explain 
both the presence and the absence of crustal fields in dif- 
ferent locations on the Moon and in rocks of different ages. 
Understanding the origins of the lunar magnetic fields has 
relevance to the understanding of terrestrial planets and their 
evolution, with implications for the magnetization of smaller 
bodies, either through dynamo generation or other means of 
magnetic field generation or enhancement. 

[ 9 ] In this study, we follow the method of Whaler and 
Langel [1996] and Wlraler and Purucker [2005] to produce 
a spatially continuous magnetization model for the lunar 
crust as a global layer with constant thickness, using inverse 
techniques based on those developed for terrestrial and 
Martian satellite data sets. As well as being able to describe 
the observed magnetic fields, magnetization models should 
be able to provide information on the properties of the 
magnetized sources, namely the magnetization contrast, 
source thickness and dimensions. 

2. Data 

[ 10 ] Crustal magnetization models require crustal mag- 
netic field measurements, which must be isolated from the 


total magnetic field measured by the MAG. The magneti- 
zation models presented here are produced using the radial, 
theta and phi components of the magnetic field synthesized 
from a degree 150 spherical hannonic model of the crustal 
field [ Purucker , 2008], and also the altitude normalized 
(to 30 km, using an equivalent source dipole technique) 
radial component of the crustal magnetic field isolated by 
Purucker [2008]. These data utilize the low altitude (<66 km, 
average 30 km), high inclination phase of the MAG data 
from the lunar wake and magnetotail times where the solar 
wind plasma densities are low and external field variations 
minimal. Altitude normalized, rather than original MAG 
data, were used because the normalization process was part 
of isolating the crustal field from remaining smaller scale 
external fields, and MAG data at the original altitudes (not 
normalized or modeled) are noisy. Description of the isola- 
tion of the crustal field from the original MAG data and the 
production of the spherical harmonic model are found in 
Purucker [2008], along with links to the data sets. 

[ 11 ] Since the production of these models, a new spherical 
harmonic model for the lunar crustal field has been produced 
by Purucker and Nicholas [2010], which supersedes the 
maps of Purucker [2008], However, the degree by degree 
correlation between the models of Purucker [2008] and 
Purucker and Nicholas [2010] is higher than 0.75 for 
degrees 15-170, and the magnetic field predictions show a 
linear correlation coefficient of 0.9. While an error was 
identified in the isolation process of Purucker [2008], 
reducing the magnitude of the radial component, this would 
not make a significant difference to the models described 
here, and the inferences from them. 

[ 12 ] Spherical harmonic terms up to degree 150 (consid- 
ered robust [Purucker, 2008]) were used to reproduce the 
vector components of the magnetic field at 30 km altitude 
(defined above the surface of a spherical Moon of radius 
1737.1 km) at evenly distributed points. Data locations were 
assigned using a ‘polar coordinate subdivision’ [Katanforoush 
and Shahshahani, 2003], placing equally spaced points on 
equally spaced (1°) latitudes, resulting in 35812 data loca- 
tions with mean spacing — 1.1°. This gives crustal field 
values up to ±18 nT, which should be considered a mini- 
mum estimate of the strength of the field at a given location, 
due to the smoothing out of smaller scale, higher magnitude 
features in the spherical harmonic representation. The over- 
all noise level of the synthesized and processed MAG data 
is dominated by the external magnetic field environment, 
and amounts to ~0.5 nT. 

[ 13 ] The radial field component of the crustal field isolated 
by Purucker [2008] retains higher amplitude unsmoothed 
information from each satellite pass in the lunar wake and 
magnetotail, and though altitude normalized, contains dis- 
continuities due to upward or downward continuation from 
the initial data altitude. A sub-set of these data was selected 
at points closest to the locations of a polar coordinate sub- 
divided grid with —1.1° mean spacing. Subsequently, points 
selected closest to the same grid shifted by 0.25° in latitude 
and longitude were used to assess the stability of the mag- 
netization solutions with respect to the input data, and pro- 
vide an upper bound on the uncertainties in the deduced 
magnetizations. 

[ 14 ] The scalar estimates of the surface field from the 
Electron Reflectometer are not used to produce these 
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models, but the (equal angle) 1° binned ER data are com- 
pared to the surface field predictions from these models, as 
an independent validation of the ability of these models to 
reproduce magnetic field data. 

3. Modeling 

[ 15 ] The inverse modeling method follows that of Whaler 
and Purucker [2005], using a minimum norm strategy to 
find a magnetization solution with the minimum root-mean 
square (RMS) magnetization amplitude for a global layer 
with constant thickness. This avoids magnetizations stronger 
than absolutely required to fit to the data. An alternative 
norm to minimize would be a measure of the spatial gradient 
of the solution [e.g., Constable et al., 1987]. However, this 
would smooth out any rapid changes in magnetization 
expected from sharp geological boundaries, and edges of 
features associated with magnetization or demagnetization 
processes. 

[ 16 ] While we do not expect a magnetized layer of uni- 
form thickness, given the lack of information about possible 
source locations and dimensions for the Moon (and the lack 
of clear correlations between magnetic anomalies with 
Clementine gravity and topography data [ Frey et al., 2002]), 
a magnetized layer with constant thickness is not an unrea- 
sonable assumption. For the Earth, the whole of the crust 
down to the Curie isotherm or the crust/mantle interface 
(whichever is shallower) can be considered to be magnetized 
in the presence of the global magnetic field. However, for 
the Moon the actual region of the crust which is magnetized 
will depend on the extent of the magnetizing fields, and the 
types of sources. 

[ 17 ] The thickness of this global layer, extending from the 
surface, was varied from 5 km to 60 km, to reflect the pos- 
sible magnetization thicknesses for the Moon; magnetization 
may exist as a thin layer of impact ejecta or mare basalt, or 
within the entire thickness of the crust, estimated from Lunar 
Prospector and Clementine gravity and topography data to 
be 49 ± 16 km on average [ Wieczorek et al., 2006], Halekas 
et al. [2002] used the demagnetization signatures of large 
(>50 km) craters to place an upper limit on the thickness of 
the lunar magnetized layer at ~50 km [ Halekas et al., 2002], 
Following Voorhies et al. [2002], the global magnetic crustal 
thickness was estimated from the spectral content of the 
lunar crustal field to be 25-30 km (see Appendix A). 
Applied to satellite data, our modeling technique is not 
expected to be able to resolve depth variations of magneti- 
zation. However, since we used data normalized to, or syn- 
thesized at, only 30 km altitude, varying the thickness of the 
magnetized layer may indicate the most appropriate average 
global thickness of the magnetized layer, which may in turn 
offer some constraints on the types of sources (e.g., volcanic 
or ejecta) or on the magnetizing field. 

[is] The continuous magnetization vector M(s), for three 
components of magnetization is expressed as a linear com- 
bination of Green’s functions G(s,r ( ), with expansion coef- 
ficients a relating the magnetization at a location s within 
a layer to N magnetic field observations at data locations r, 
[ Parker et al., 1987; Share et al., 1982]: 

N 

M(s) = ^>,G(s,r,) (1) 

i=l 


The Green’s function relating magnetization to the />th 
component of the observed magnetic field is given by 


with t| /,; the unit vector in the direction of the local orthog- 
onal components. Application of the gradient operators with 
respect to positions r and s give basis functions (and hence 
magnetization) that decrease with depth into the magnetized 
layer, inversely proportional to |r, — s| 3 . 

[ 19 ] The set of magnetic field observations B, of single or 
multiple components, can be represented as a system of 
linear equations: 

B = Ea (3) 

where T is the ‘Gram matrix’ containing the volume inte- 
grals of all the possible pairs of dot products of the basis 
functions for the data points ordered in vector B, and a is 
the vector of coefficients of the basis functions given in (1). 
A solution for a with minimum magnetization amplitude 
for a given fit to the data is obtained from [Shure et al, 
1982] 

B = (r + AI)c* (4) 

where I is the identity matrix and A, known as a damping 
parameter, controls the relative importance of the fit to the 
data versus the amplitude of the magnetization solution. 
As the damping parameter is increased, the magnitude of 
magnetization, and any features of the model relating to 
noise in the input data, are expected to be reduced at the 
expense of the fit to the data. 

[ 20 ] The system of equations (4) is too large to be solved 
directly for the data sets considered here. For a magnetiza- 
tion distribution confined to a layer of constant thickness, 
the forms of the elements of the Gram matrix, given in 
Jackson [1990] and Whaler and Langel [1996], can be 
expressed as functions of the cosine of the angular separa- 
tion of data points. For large angular separations, they 
become very small and the Gram matrix can be made sparse 
by treating as nonzero only matrix elements with absolute 
values above a stated threshold. Following the method of 
Whaler and Purucker [2005], a solution for the expansion 
coefficients can be found using an iterative conjugate gra- 
dient method for sparse matrix systems [ Purucker et al., 
1996], with Jacobi scaling. The solution was deemed to 
have converged when equation (4) was satisfied to a toler- 
ance of <10- 10 , within 500 iterations. Failure to converge 
within 500 iterations would indicate an ill-conditioned sys- 
tem, requiring a larger damping parameter. 

[ 21 ] Alongside the solution vector a, the sum of squares 
of residuals (SSR) was calculated from [Shure et al., 1982] 

SSR = \ 2 a T a (5) 

(more usefully expressed as the RMS misfit), and the solu- 
tion norm ||M||, representing the RMS magnetization 
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Table 1. Model Outputs for a 30 km Thick Magnetized Layer d 


Damp. A 

SSR (nT 2 ) 

Norm (A/m) 

Misfit (nT) 

Iterations 

Br 

Correlation Coefficients 
Bt 

Bp 

0 

0 

2228 

0 

30 

1 

1 

1 

100 

0.300 

2223 

0.0017 

30 

1 

1 

1 

1000 

27.44 

2175 

0.0160 

26 

1 

>0.999 

>0.999 

5000 

500.2 

2010 

0.0682 

22 

0.999 

0.999 

0.998 

10000 

1536 

1858 

0.1196 

18 

0.998 

0.997 

0.996 

50000 

13090 

1246 

0.3490 

11 

0.985 

0.983 

0.98 

100000 

25409 

909.4 

0.4863 

8 

0.976 

0.974 

0.969 


“Models produced from three components of the magnetic field synthesized from the spherical harmonic model of Purucker [2008], Linear correlation 
coefficients [e.g., Press et al., 1992] are between the predicted and input data. 


amplitude for the entire volume of the magnetized layer, was 
also calculated: 


|M|| 


N 


J M ■ MdV = t/cdTa 

V 


( 6 ) 


The three components of the continuous distribution of 
magnetization obtained by substituting the solution for a 
from (4) into (1) were evaluated on grids with 0.5° and 0.25° 
spacing, at a specified depth within the magnetized layer. To 
allow a more meaningful comparison of models with dif- 
ferent layer thicknesses, we calculated the vertically inte- 
grated magnetization (VIM, in units kA). 

[ 22 ] Ideally the crustal magnetization models would be 
able to adequately account for both the MAG data at altitude 
and the ER data at the surface. The magnetization solutions 
were also used to predict the magnetic field at the surface 
locations of the 1° binned ER data for comparison. How- 
ever, the surface fields predicted from the magnetization 
model are not expected to agree in magnitude with the ER 
data due to the differences in the sensitivity of the MAG and 
the ER technique. The ability of the magnetization model to 
predict the ER data acts as a test for the presence of anni- 
hilators of the magnetic inverse problem [Maus and Haak, 
2003], a source of non-uniqueness [ Runcorn , 1975]. Mag- 
netic annihilators contribute to the solution without produc- 
ing magnetic fields at the observation locations; however, 
they may produce nonzero magnetic fields at other locations. 


4. Results 

[ 23 ] Solutions for models from three components of the 
synthesized magnetic field with a range of layer thicknesses, 
and damping parameters between 0 and 100,000, converged 
within 80 iterations or fewer. Fewer iterations were required 
for models with larger damping parameters. Models from the 
radial field component data required higher damping para- 
meters than those used for the models from the synthesized 
data. 

[ 24 ] The robustness of the solution to the number of non- 
zero elements in the sparse Gram matrix was tested, and for 
a 30 km thick layer, threshold values from 0.1 to 10 retained 
from 9.18% to 0.988% of the Gram matrix elements 
respectively with very little difference (to within less than 
0.05%) in the SSR and the solution norm for each of these 
models. A threshold value of 1.5 was used for the models 
with layer thickness from 8 km to 60 km, retaining between 


1 .3% and 3.7% of the Gram matrix elements. Thinner layers 
needed lower threshold limits to converge and are deemed 
less reliable due to potential instabilities in the computations. 

[ 25 ] Model solutions from three components of the syn- 
thesized magnetic field fit the input data to within the esti- 
mated noise level (see Table 1 for a 30 km thick layer), with 
globally averaged misfits of <0.5 nT, and high global linear 
correlation coefficients (>0.94) between the predicted and 
input synthesized fields. Solutions from the radial field 
component data resulted in higher misfits due to the greater 
range in the input data values, so while mentioned here are 
not discussed in as much detail as the models from the 
synthesized data. 

[ 26 ] While a convergent solution was found with zero 
misfit when no damping was used with the three compo- 
nents of the synthesized field, such a model is likely to 
contain features which relate to noise in the data. Damping 
parameters between ~10,000 and 50,000 (comparable to the 
magnitudes of the diagonals of the Gram matrix) give 
models with similar fits to the data having lower RMS 
magnetizations, so represent models which optimally mini- 
mize both the residuals and the magnetization. 

[ 27 ] Discussions of the possible thickness of the magne- 
tized layer are based on model solutions from the three 
components of synthesized data, for layers with thicknesses 
of 10 km, 30 km and 50 km with the same misfit (of 
~0.2 nT). The component and total magnetization solutions 
for these models, calculated at the surface, can be found in 
the auxiliary material. 1 These models have the same number 
of iterations despite their different damping parameters and 
thicknesses, and give very similar, and high, correlation 
coefficients and similar VIM (Table 2). The total magneti- 
zations for these models at the top and bottom of the layer 
are shown in Figure 1, showing that the pattern (and 
direction) of magnetization remains constant with depth but 
decreases in amplitude. The 10 km thick layer requires the 
lowest damping parameter to give the same fit to the 
data as the 30 km and 50 km thick layers, and shows a 
higher solution norm for the same damping parameter. The 
‘intermediate thickness’ model, of 30 km, has the smallest 
misfit and smallest RMS magnetizations for most damping 
parameters: models with thicker and thinner layers require a 
larger RMS magnetization to give the same misfit as that for 
the 30 km thick layer model. 


'Auxiliary materials are available at ftp://ftp.agu.org/apend/je/ 
201 lje003944. 
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Table 2. Model Outputs, Including Fit to the ER Data, for Selected Models With Different Thickness Magnetized Layers 3 


Correlation Coefficients 


Thick, (km) 

Damp. A 

SSR 

(nT 2 ) 

Nonn 

(A/m) 

Misfit 

(nT) 

ER Misfit 
(nT) 

VIM 

(kA) 

Br 

Bt 

Bp 

10 

15000 

4581 

1991 

0.207 

16.5 

0.094 

0.997 

0.997 

0.996 

30 

20000 

4227 

1638 

0.198 

17.0 

0.139 

0.994 

0.993 

0.991 

50 

16000 

4164 

1761 

0.197 

17.2 

0.208 

0.993 

0.989 

0.986 


“Models produced from the three components of the magnetic field synthesized from the spherical hannonic model of Purucker [2008], Linear 
correlation coefficients [e.g., Press et al., 1992] are between the predicted and input data. 


[ 28 ] For magnetization confined to layers with thicknesses 
from 10 to 50 km, magnetization amplitudes of at least 
75 mA/m are required to produce the strongest observed 
magnetic fields. The average surface magnetization is very 
weak at 5.4 mA/m, while the crustal average is 4.9 mA/m 
(for the 30 km thick layer), reflecting a modest reduction in 
magnetization amplitude with depth. Average values in 
magnetized regions are ^30 — 40 mA/m, while maximum 
values in the most strongly magnetized areas such as the 
Reiner Gamma region are up to 120 mA/m. 

[ 29 ] Figure 2 shows the three component magnetizations 
calculated at the surface for a model from the synthesized 
data with a 30 km thick magnetized layer, selected to have 
misfit 0.2 nT. The magnetizations in the radial direction are 
the strongest, and the magnetizations in the theta direction 
are the weakest. The magnetizations vary in direction over 
small regions, mirroring the directions and relative ampli- 
tudes of the vector components of the crustal field. The 
magnetization features seen in this model are present in 
those for all thicknesses and damping parameters, but with 
different amplitudes. 


[ 30 ] Magnetization solutions from the radial field compo- 
nent on the global scale are almost identical to the solutions 
obtained from the synthesized data in Figures 1 and 2, and 
are hence not shown separately. However, on a smaller 
scale, the solutions from the radial data show finer scale 
features at the edges of magnetized regions. Magnetization 
distributions for models from the radial data on shifted grids, 
while having the same solution nonn and SSR (to within 
1%), showed visible differences to each other in the finer 
scale detail of the patterns of magnetization, highlighting an 
instability with respect to the inversion of unsmoothed data. 

[ 31 ] Predictions of the magnetic field observations (using 
equation (3)) from both the models from the synthetic data 
and the radial data show that the pattern of most of the 
small scale and strong features of the field are reproduced, 
but with reduced magnitudes. The absolute residuals, defined 
as \data\ — \predictions\, shown in Figure 3 for the model 
from the synthetic data, show this under-prediction of the 
data as positive residuals, and also over-prediction (negative 
residuals) around the edges of the strongest anomalies. The 
largest residuals, of up to ±5 nT (higher for the more noisy 



0 


5 


15 30 50 75 


Magnetization (mA/m) 

Figure 1. The total magnetizations at (top) the surface and (bottom) base of the magnetized layer 
for models with magnetized layers of (a) 10 km thickness with A = 15,000, (b) 30 km thickness with 
A = 20,000, and (c) 50 km thickness with A = 16,000, from three components of the magnetic field syn- 
thesized from the spherical hannonic model of Purucker [2008] (Lambert projection centered on the far 
side, 0.5° gridded points). 
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Magnetization (mA/m) 


Figure 2. Three component magnetizations at the surface 
for the global magnetization model (with a 30 km thick layer 
and A = 20,000) from three components of the magnetic field 
synthesized from the spherical harmonic model of Purucker 
[2008] (Lambert projection centered on the far side, 0.5° 
gridded points). 


radial data), occur where the measured field is maximum, 
and for large damping parameters. 

[ 32 ] Power spectra of the radial component of the mag- 
netic field predictions (Figure 4) are weaker than the power 
of the input field, showing the reduced magnitudes. The 
resulting spectra for predictions from the selected model 
with damping parameter 20,000, and an ‘over-damped’ 
model (A = 50,000) show the good correlation between the 
predictions from these magnetization models and the input 
data up to about harmonic degree 100. The decrease in 
power at higher harmonic degrees results from the use of the 
damping parameter in the inversion. A similar fall off in 
power for predictions from highly damped models is seen 
for the predictions from the models with the radial data, 
where higher damping parameters were used. 

[ 33 ] The preferred models with 10 km, 30 km and 50 km 
thick magnetized layers have very similar global distributions 


of residuals. However, there are more small residuals in the 
model with 10 km thickness, and larger residuals in the 
model with 50 km thickness, despite these models having 
the same total misfit. The 10 km thick model fits some of 
the strong anomalies better than the other thicknesses, 
though other regions showed very little variation in residual 
distribution with model thickness. 

[ 34 ] When used to predict the surface magnetic fields, the 
models from the radial component data, and the 10 km thick 
magnetized layer from the synthesized data, gave slightly 
better fits to the ER data (Figure 5 for 10 km thick layer). 
The predictions and the ER data agree well spatially, with 
models of all magnetized layer thicknesses giving the same 
correlation coefficient with the ER data (to 3 significant 
figures) of 0.639. However, the actual agreement of the 
ER and predictions varies for each region, with the stron- 
gest total surface fields being under-predicted, and the 
weaker field region values being over-predicted. The largest 



Residual (nT) 


Figure 3. Absolute residuals for the three components of 
magnetic field from the selected global magnetization 
model with 30 km thickness. Positive residuals signify 
under-prediction, negative residuals signify over-prediction 
of the input data. (Lambert projection centered on the far 
side, 1° gridded points). 
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Figure 4. Comparison of the spatial magnetic power spectra, Rn (at the surface), of the synthetic radial 
data (blue crosses) [ Purucker , 2008] and radial predictions from the selected model with A = 20,000 (green 
circles) and an ‘over-damped’ model (A = 50,000, red diamonds). 


differences between the data and predictions occur in the 
regions of very strong (>10 s nT) or very weak (<few nT) 
surface magnetic field. 

5. Discussion 

[ 35 ] The magnetization distribution and amplitude 
remains considerably robust for a range of damping para- 
meters and thicknesses. Magnitudes differ only by a factor 
of <4, allowing constrained estimates of likely lunar crustal 
magnetization, unlike the spatially continuous models for 
the Martian crust which saw the RMS magnetization vary 
by a factor of 20 for models with misfits different by only a 
factor of 2 [Whaler and Pumcker, 2005]. An upper value of 
the intrinsic model uncertainty in the magnetization values 
in our models can be estimated from the differences 
between the magnetization solutions produced when the 
radial input data grid was shifted. For the Reiner Gamma 
region, the maximum differences are 20-25 mA/m for the 
largest values of magnetization. While this can be more than 
100% of the magnetizations at some locations, such large 
differences result from the slight shifting of the resulting 
magnetization distribution with each input data set, and the 
mean difference of ~4 mA/m, for regions with magnetization 
>5 mA/m, is a better estimate of the model uncertainty. 

[36] There is ‘one to one’ correspondence of the directions 
of magnetization and the vector components of the crustal 


field, despite the inversion not fixing source locations, 
meaning that the magnetization directions in these models 
are unlikely to be representative of true lunar magnetization 
directions. The horizontal extent of the sources, though 
coincident with the magnetic field feature, may be more 
compact than these magnetization models show due to the 
horizontal spreading of the magnetization solution across the 
modeled layer. As a result, there is no additional information 
about signatures associated with particular lunar anomalies. 

[ 37 ] An inversion technique minimizing spatial variation 
rather than amplitude may have resulted in a solution with 
less directional variation. The magnetization modeling 
method used here however, is not conducive to producing 
large regions with uniform magnetization direction, so cal- 
culation of the magnetization inclination and declination is 
not meaningful. In addition, care is required in the inter- 
pretation of the magnetized sources in some regions where 
linear looking features in one or more magnetization com- 
ponent may not be linear when the total magnetization, or 
another magnetization component is considered. Such cor- 
respondence between the components of magnetization and 
the components of the field is also observed in crustal 
magnetization models at Mars and the Earth [Whaler and 
Purucker, 2005, 2010]. 

[38] While magnetization directions from these models 
are to be treated with caution, the magnetization strengths 
appear to be reliable. The magnetization strengths of these 
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Figure 5. Comparison of (left) the mid latitude 1° binned ER data with (right) the surface total field pre- 
dictions at the same locations from the global model with a 10 km thick layer and A = 15,000. 
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models satisfy the analytical bound to account for the 
observed magnetic fields provided from ideal body theory 
[ Parker , 2003] and are consistent with measured and esti- 
mated values of lunar magnetization. Magnetization values 
of 0.5-200 mA/m, as mass normalized dipole moments 
(for magnetizations on a grid with 0.25° spacing) are 
~4 x 10” 8 — 1 x 10~ 4 Anrkg -1 , comparable to the 
remanent magnetizations measured for the Apollo samples at 
10 -7 — 10 _1 Am 2 kg _1 [ Fuller and Cisowski, 1987], and 
differentiated meteorites (~10“ 3 Anrkg -1 ) [Kivelson et al, 
1993]. Mitchell et al. [2008] estimated the magnetization ofa 
uniformly magnetized impact melt in the Imbrium basin 
based on the weak fringing fields of craters within this basin 
at 10~ 6 Am 2 kg _1 , and this model shows similar magnitudes 
in that region. 

[ 39 ] The minimization of the solution norm means that 
these magnetization magnitudes are minimum estimates of 
lunar crustal magnetization. Like the magnitudes of the 
crustal fields, the lunar crustal magnetizations are very much 
weaker than those on Earth and Mars for layers of compa- 
rable thickness. In fact, the highest magnetizations on the 
Moon of ~0.2 A/m are barely significant compared to typical 
magnetizations observed on Earth of 1-5 A/m [Wasilewski 
and May hew, 1982]. While thinner layers and smaller scale 
features are likely to have much stronger magnetizations 
than the maximum values obtained here for a layer with 
constant global thickness, lunar magnetizations are unlikely 
to be comparable in magnitude to those on Earth and Mars 
due to the differences in magnetic mineralogy, the strength of 
the magnetizing fields, and the processes of magnetization. 

[ 40 ] The global misfits for these models are within the 
noise levels of the input data, and the difference in misfits 
for models with magnetized layers in the range 10-50 km 
thick are small; thus there is no clear preferred layer thick- 
ness. However, the smaller misfits and RMS magnetizations 
for the model with a 30 km thick layer, and the persistence of 
magnetization features at depths greater than 10 km but less 
than 50 km (Figure 1), suggests that globally the magnetized 
sources of the lunar crust can be represented by a layer 
~30 km thick. It is also possible to estimate a magnetic 
layer thickness using Parker's [2003] ideal body method in 
a reverse sense; if the magnetizations are the minimum 
required to give the fields at altitude, Parker's [2003] 
expressions show that layers could be as thin as 3 km for 
the strongest magnetizations, but the average magnetiza- 
tions give thicknesses of 10-15 km. It is not expected that 
all the magnetized regions will have the same thickness, and 
there are likely to be regions over the globe which can be 
better described by thinner or thicker layers of magnetized 
crust (for example the thicker crust of the lunar farside). 
Indeed, the residuals over some regions, and the agreement 
of predictions with surface field estimates, are improved 
slightly if a thinner magnetized layer (~10 km) is used. 

[ 41 ] The disagreement with some of the surface data 
results from the difference in sensitivity of the MAG and 
the ER technique to sources with different wavelengths, and 
the use of the damping parameter in these inversions. The 
strongest surface features are under-predicted by the model, 
as the ER technique is more sensitive to smaller wavelength 
sources than the MAG. Over-prediction of the weaker sur- 
face fields is also observed with the spherical harmonic 


model, and may occur as a result of the amplification of 
noise in the MAG data when effectively downward con- 
tinued to the surface, or demonstrate the variation in the 
response of the ER technique to the very small wavelength 
features anticipated from the Apollo surface data. 

6. Conclusions 

[ 42 ] Spatially continuous global crustal magnetization 
models for the lunar crust have been produced using the 
method of Whaler and Purucker [2005]. These are the first 
global models of lunar crustal magnetization, and have 
magnetization solutions with minimum RMS magnetization 
amplitude for a given fit to the data, controlled by a damping 
parameter. Models were produced from both the radial 
component of the MAG data from the wake and tail regions, 
and from the three components of the magnetic field syn- 
thesized from a degree 150 spherical harmonic model of the 
lunar crustal field. The resulting distributions of crustal 
magnetization produced from both data sets show similar 
patterns, demonstrating robustness. The magnetization 
models can account for many of the features of the observed 
magnetic fields at altitude, and for some of the surface fields 
as measured with the ER, though amplitudes are often not 
reproducible. 

[ 43 ] Surface predictions from these magnetization models 
provide a better prediction of some of the weaker surface 
fields than the degree 150 spherical hannonic model. How- 
ever, as a single layer thickness globally will not be appro- 
priate everywhere, and a single choice of model will fit 
individual data points differently, predictions from this 
magnetization model fail to capture some of the field infor- 
mation, particularly on small spatial scales (higher harmonic 
degrees), contained in the unmodeled data or spherical har- 
monic model. 

[ 44 ] The suggested 30 km thick magnetized layer is con- 
sistent with the spectral thickness calculated, and represents 
a global compromise of the fit to the data and RMS mag- 
netization. For example, thinner layers provided a better fit 
in the Reiner Gamma region, compatible with theories that 
the magnetized source in this region is a thin near-surface 
layer of impact ejecta [Hood et al., 2001; Nicholas et al., 
2007] (though the layer would likely be thinner than the 
thinnest modeled here). Models with layer thicknesses 
>10 km require magnetization magnitudes that are compa- 
rable to those measured in the Apollo samples, demon- 
strating that minimum magnetization models from satellite 
data provide useful bounds on the magnetization of the 
lunar crust, and suggest that the magnetic mineralogies of 
the Apollo samples are representative of the lunar crust as a 
whole. 

[ 45 ] The magnetization strengths suggest that magnetizing 
fields from a few /iT (as deduced for some Apollo samples) 
to more than 40 /iT (for shock magnetization [Gattacceca 
et al., 2010]) might have been present on the Moon. The 
persistence of these magnetization strengths at depth (>10 km) 
in some regions requires the magnetizing field to exist in, or 
be able to penetrate into, the lunar crust, in addition to the 
mechanism for producing the remanence. As transient fields 
generated in impacts are expected to penetrate into the lunar 
crust [ Hood and Artemieva, 2008], magnetization at depth 
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does not distinguish impact generated fields from a dynamo 
field unless magnetization was acquired during cooling, 
requiring a field steady over long periods of time. However, 
in the absence of field amplification, a dynamo field may not 
be strong enough to explain some of the stronger magneti- 
zations, and thus a combination of different fields and 
mechanisms may be needed to explain the lunar crustal 
fields. 

[46] While this modeling method did not permit magneti- 
zation directions to be determined, some of the spatially 
rapid directional variation in these magnetization models 
may be real features of lunar crustal magnetization, sug- 
gesting that the directions of the magnetizing fields may 
have varied on small spatial scales. Such small scale varia- 
tion within a magnetizing field is more likely to be associ- 
ated with impact generated magnetic fields than a dynamo 
field. 

[ 47 ] With the excellent coverage of the LP (and now also 
Kaguya) MAG data [e.g., Tsunakawa et al., 2010], and new 
models of the internal magnetic field [Purucker and 
Nicholas, 2010], it is not necessarily the resolution or cov- 
erage of the crustal field measurements which are limiting 
the conclusions on the sources of lunar crustal magnetism 
(though surface surveys such as that from a lunar rover 
would of course provide valuable information); rather, as in 
the terrestrial case also, it is the deciphering of the exact 
sources of the fields which is still lacking. 

[48] Ultimately, a better understanding of the magnetiza- 
tion of the lunar crust and the sources themselves will come 
from regional magnetization models, more lunar samples, 
and new geophysical data. Ground truth magnetization 
directions and rock compositions for regions like Reiner 
Gamma and impact antipodes, coupled with, for example, 
seismic and/or radar surveys, may constrain the true sources 
of these anomalies and their magnetizing fields. Re-analysis 
of the magnetic properties of the Apollo samples, and 
experimental studies of shock remanent magnetization [e.g., 
Gattacceca et al., 2007, 2008, 2010] is ongoing to improve 
the understanding of the magnetic and paleomagnetic prop- 
erties of the lunar rocks. 

Appendix A 

[ 49 ] Voorhies et al. [2002] used the theoretical power 
spectrum expected from a large number of uncorrelated 
dipoles, with random magnetic moments, scattered on a 
sphere within a spherical body to find the best fitting source 
radius for different sections of the observed power spectra 
for the crustal field of Earth and Mars. Applying this tech- 
nique to the Moon, the best fitting power spectrum for 
degrees 1 to 150 is shown in Figure Al with the observed 
power spectrum. The higher degree terms fit the theoretical 
model better than the lower degree terms, which are perhaps 
dominated by noise from unmodeled external fields. 

[ 50 ] The best fitting source sphere radius of 1711.4 km 
gives a globally averaged source depth of 25.7 km. The 
globally averaged source depth varies from ~14 to 30 km 
depending on which parts of the spectrum are fit. The esti- 
mated source depth from this technique is considered to be 
closer to the layer thickness than the layer half thickness as 
theoretical spectra may overestimate the true depth of sour- 
ces that are correlated over scales larger than the dipole 



Figure Al. The spatial magnetic power spectra, Rn, for the 
lunar crustal field up to degree 150 at the surface (blue 
crosses) with a theoretical curve fitted through it (red line) 
expected from a spherical shell of randomly oriented dipoles. 


spacing [Voorhies, 2008]. Although this ‘spectral thickness’ 
does not take into account any correlation of the sources, or 
variations in layer thickness which would be expected for 
some regions of the Moon, a layer thickness of 25-30 km 
serves as a first estimate for the global average magnetic 
crustal thickness of the Moon. 
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